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A novel method for the precise identification and determination of the energies that contribute in 
the spectral decomposition of lattice correlators is presented. The method is based on statistical 
concepts and it rehes heavily on simulation techniques. The r\c correlator is analyzed within this 
method and the results obtained are compared to a previous analysis based on Bayesian statistics. 
An analysis of the nucleon local two-point correlators leads to the identification of the excited 
states in the positive and negative channels. A discussion on the Roper is included. 
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1. Introduction 

A significant component of hadronic pliysics researcli focuses on the understanding of the ex- 
citation spectrum of hadrons and in particular that of the proton [[T]]. On the experimental side, 
high quality data are becoming available largely due to the advances in accelerator and instrumen- 
tation technologies. Lattice QCD simulations using new algorithms and faster computers yield 
high-precision results close to the chiral limit, where chiral perturbation theory can be reliably ap- 
plied thus yielding physical results that can be compared to experiment. However, achieving a high 
precision necessitates a robust analysis of simulation data keeping systematic errors under control. 
Such an analysis typically involves a fitting procedure. The quantities of interest are masses of 
hadrons, decay constants, form factors and other hadronic matrix elements. In most cases one uses 
a subset of the lattice data and truncates the theory to a small number of fitting parameters. For 
example, in the case of two-point functions one discards simulation results involving short times. 
Such a truncation eliminates important information about excited states encoded in these two-point 
functions. The goal of the current study is to apply a method to extract the maximum information 
from lattice measurements. 

Various approaches have been proposed to ameliorate the problems of the conventional fitting 
procedure most of which are based on the Maximum Entropy Method ^ Q ^. In this work we 
present a new method that relies solely on ;^^-minimization with an unbiased evaluation of errors. 
It has been developed and applied in the context of the analysis of electroproduction data in the 
nucleon resonance region [^] but, as demonstrated here, it is of general applicability. The method 
has as a minimal requirement that the parameters to be determined are linked in an explicit way to 
the data. There is no requirement that this set of parameters provide an orthogonal basis. Moreover 
these parameters can be subjected to explicit constraints. The main advantages of the method is that 
it requires no prior knowledge (priors) other than the spectral decomposition of lattice correlators 
nor any assumption on the level of truncation. The data set alone determines the information that 
can be extracted. 

In order to demonstrate the method we analyze the T7c-correlator to extract the three lowest 
states [0] and compare the results of our method to the values extracted using the method of Ref. [Q] 
on the same data. Local nucleon correlators on standard lattices are easy to produce and are readily 
available. We show that the method is well suited to extract the excited states of the nucleon in 
both the positive and negative parity channels. Using dynamical twisted mass fermions as well as 
dynamical Wilson fermions and two different nucleon interpolating fields we discuss the positive 
parity excitation of the nucleon and its connection to the Roper as a function of the pion mass. 

2. The method 

The method, referred to as "AMIAS" [^, relies only on the ergodic hypothesis, namely that 
any parameter of the theory can take any possible value allowed by the theory and its underlying 
assumptions. The probability of this value representing reality is solely determined by the data. 
We assume that all possible values are acceptable solutions, but with varying probability of being 
true. If n parameters are needed to describe the data, then each j-set is a solution. We 

assign to each such set a ;^^-value and a probability. We then construct an ensemble of solutions 
and assume that this contains all solutions with finite probability. The probability distribution for 
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any parameter assuming a given value is then the solution. An ensemble of solutions can be defined 
as the collection of solutions, which are characterized by ^ Xmin + ^i'^^ 'tf usually taken to 
be a constant equal to the effective degrees of freedom of the problem. For a solution set j we 
compute X^if^j j) = Lili L/{ ^^''^'''^wk{t')^^^^^ J"^' ''^here Vk{ti) and Wk{ti) are M sets of measurements 
and errors respectively. For lattice applications this set of measurements is defined on temporal 
lattice slices t,- and f{tj,{B}^) is the spectral decomposition for the relevant correlator written in 
terms of the n parameters {B}j. In what follows we further develop the methodology through the 
specific problem of extracting the excited states from two-point correlators. The Euclidean time 
correlator C{t) of an interpolating operator J{x,t) and its spectral decomposition for zero three- 
momentum is 

oo 

Cit)=Y^<Jix,t)j\0,0)>=Y^Aie-""' , (2.1) 

X /=0 

and therefore for each j-set f{t,-,{A,m}) = Lf=o^'^ '"'''> where L is the highest excited state that 
the data are sensitive on and it is determined by the method. The exponential dependence is cor- 
rect if we neglect boundary conditions (b.c). For antiperiodic b.c. a meson correlator is sym- 
metric about the mid point of the temporal lattice extension and each exponential is modified to 
coshm(f — r/2). We start by assuming a maximum number of L excited states in the spectral de- 
composition of the correlator and select a suitable range of values for each of the parameters A/ > 
and tfiQ < mi < m2 < ■ ■ ■■ These ranges define the "phase volume" which the Monte Carlo method 
explores. Within the chosen range, we uniformly select a value for each parameter. Considering 
L excited states, i.e. L + 1 exponentials with n = 2 + 2L parameters, we evaluate the x^ (2 + 2L, j) 
corresponding to the particular solution using our lattice data. We repeat this procedure a large 
number of times, typically a few hundred thousand, generating an ensemble of solutions. The 
method does not depend on the choice of L, provided that it is sufficiently large. In practice it is 
chosen so that the derived results do not change if instead L + 1 is used as a maximum cutoff. The 
computational time required depends critically on the choice for the phase volume to be explored. 
For the implementation described here it scales like a power in the number of parameters. However 
this can be reduced to about linear if a Markov chain is used to generate the ensemble. It is obvious 
that the overwhelming majority of the solutions generated in this fashion are characterized by very 
large ;(^^-values. However, as shown in Fig. [l] a saturation of "good values" is achieved as the phase 
volume is enlarged that remains unaltered by exploring a wider range of values for the parameters. 

The behavior of a given parameter A/ or m/ in the ensemble of solutions can be visualized by 
plotting the x^ versus the value of the parameter. In Fig. |2| we show such a plot for mo, the ground 
state mass of a two-point correlator. The sensitivity of the data on this parameter becomes explicit 
and quantifiable by applying successive cuts on the x^- values and constructing out of the selected 
population of solutions histograms, as shown in Fig. ^. For this parameter on which the data are 
sensitive on, the distributions have a well defined maximum independent of the cut whereas, as 
expected, the width depends on the ;^^-cut applied. If the data show no sensitivity on a given pa- 
rameter then the ensemble of solutions for this particular parameter yields a uniform distribution 
independent of the ;^^-cut.Therefore the method does not treat differently "sensitive" from "insensi- 
tive" exponential terms; they naturally emerge as such. The widths of the histograms, ranging from 
very narrow to very wide or infinite, naturally select and order the various exponentials according 
to their sensitivity to the data set. 
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Figure 1: The ;t;^-distribution as Figure 2: The ;t:^-distributions for Figure 3: The distribution of mo 
the range of values of the parame- mo of a two-point coiTelator for of a two-point correlator for dif- 
ter set {B} is increased. different x^-cnts. ferent values of Xmin- 
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Figure 4: Left: Correlations between two sensitive parameters; Right: a sensitive and an insensitive param- 
eter. 

A central issue that is properly treated in our method, is the handling of correlations. All 
possible correlations are accounted for by allowing all fit parameters to randomly vary and to yield 
solutions with all allowed values including the "insensitive" exponential terms. The visualization 
of at least the dominant correlations is accomplished in a two-dimensional scatter plot in which the 
ensemble of solutions is projected on the plane defined by the values of the parameters and color 
coded according to the ;^^-value. In Fig. Q the strong correlations between mq and mi, two of the 
parameters on which a two-point correlator depends on, are displayed as compared to the lack of 
correlations between mo and m2 on the latter of which the data do not depend on. 

3. Analysis of rjc correlator 

In this section we apply our method to re-analyze the T]c-correlator computed on seventy un- 
correlated gauge configurations on a lattice of temporal extent Nt =96 [^. These data were 
analyzed using the method of Ref. [Q] to extract the masses of the ground state and the first and 
second excited states [^. We apply the "AMIAS" method using four exponentials i.e. L = 3. In 
Fig. ^ we show the distributions for the four amplitudes and four masses. As can be seen, the three 
lowest states are determined from the data whereas the the fourth is undetermined. The values 
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Figure 5: The distribution of the amplitudes and masses for rjc- 



Figure 6: The distribution result- 
ing from the "AMIAS" ensemble 
of solutions yields the mass spec- 
trum of the rjc in lattice units. The 
widths characterize the uncertainty 
with which the method determines 
the parameters. 



that we find, starting the fit from the fourth time slice from the source, are: 1.3171(13), 1.608(9) 
and 2.010(11) in lattice units, where the errors are computed using a jackknife procedure. The 
"AMIAS" values are in agreement with the ones extracted in Ref. using the method of Ref. [Q], 
namely 1.3169(1), 1.62(2) and 1.98(22), but with improved accuracy for the excited states. Note 
that, since the mean values are determined by Xmin and the error using a jackknife procedure, the 
results are independent of the parametrization used for the function /(?,, {B}^). Compatible errors 
result from the width of the "AMIAS" distributions ||^ as can be seen in Figs. || and ^. The latter 
figure shows the Tjc-spectrum for the four states that we fitted.There is a considerable overlap of the 
excited states and this explains the strong correlations among them. It is crucial for the determina- 
tion of the mean values and the errors to increase the number of exponentials beyond the states one 
is determining. 



4. Nucleon 

As a first new application of our method we examine the excited states of the nucleon using 
local correlators that are produced easily in lattice simulations. We use two interpolating fields: 

JN{x)=e"'''{uaCY54)uc , JUx)=e"'''{ulCdh)Y5Uc . (4.1) 

We compute the two-point correlators with Jn using two degenerate flavors of dynamical twisted 
mass fermions (TMF) as well as using two degenerate flavors of dynamical Wilson fermions |^]. 
We apply our method taking L = 2 and the resulting spectrum for the positive and negative parity 
channels is shown in Fig. ^ The two lowest states can be clearly identified in both channels. In 
addition, in the case of Wilson fermions, we compute the two-point correlators using Jj^. It has 



been conjectured [10| that J'p^ has a large overlap with the Roper whereas its overlap with is 
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Figure 7: Nucleon mass spectrum in lattice units for the positive (left) and negative (right) parity channels 
using local correlators with TMF at pion mass 484 MeV on a lattice of spatial length 2.1 fm at )3 = 3.9. 




A2 M2 A2 M2 

Figure 8: Probability distributions for the amplitudes and masses in lattice units extracted from local corre- 
lators computed using Np —2 Wilson fermions at pion mass 500 MeV on a lattice of spatial length 1.8 fm 
at j3 = 6.0. Left using J^; Right using 7^. 



small. We compare the distribution of amplitudes and masses in the positive parity channel in 
Fig. H from correlators using these two interpolating fields for k = 0.1580 (mji ~ 500 MeV). One 
clearly identifies the first excited state from these low quality data. In addition, we observe that 
the state of lowest mass that is present in the mass spectrum of the correlator computed with Jf^ is 
absent when using 7^. Instead the correlator with Jj^ has a lowest state that does not show up when 
using y^v- The conjecture is that this state is the Roper. 

In Fig. ^we show the mass of the two lowest states in the positive parity channel and the mass 
of the lowest state in the negative parity channel as we vary the pion mass. Results within these two 
lattice formulations are compatible. We extrapolate TMF results on the ground state in the positive 
parity channel to the physical point using lowest order heavy baryon chiral perturbation theory, 
whereas we use a linear extrapolation in for the negative parity and first excited positive parity 
states. We show with the solid lines the best fit and with the dashed lines the error band. At the 
physical point, within the estimated errors, the ground state masses of both positive and negative 
parity channels are consistent with experiment, shown with the asterisks. 
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The mass of the first excited positive parity state 
extracted from the correlator with Jn, is too high to 
be identified with the Roper. However, the mass of 
the lowest state extracted from the correlator with 
J'j^ is in the right mass range. It is very close to the 
mass of the negative parity state and a linear ex- 
trapolation yields a value consistent with the mass 
of the Pii state, albeit with large statistical error. 
This corroborates the conjecture that this state is 
the Roper. In contrast to their different behavior 
for the ground state, both Ji\i and J'j^ interpolating 
fields yield the same mass for the first excited state 

.00 .05 .10 .15 .20 .25 .30 

in the positive parity channel as well as for the two 
m ' (GeV^) , ^ . , . . , 

" lowest states in the negative parity channels. 

Figure 9: Mass in the +ve and -ve parity channels 

of the nucleon as function of m^. 

5. Conclusions 

A novel method for identifying and extracting parameters from lattice simulation data is pre- 
sented. Applying the method to the Tjc-correlator, we reproduce the results of a previous analysis 
where the use of priors is found to be important for the identification of the two excited states. 
Our method, with no prior input, clearly identifies the two excited states with improved accuracy. 
Applying the method to the local nucleon correlators we determine the mass of the ground and the 
first excited states in the positive and negative parity channels. Our results are in agreement with 
more evolved mass correlation matrix analyzes yielding a positive parity state that can be identified 
with the Roper. 
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■ Jj^ +ve parity, Np=2 Wilson 
□ J„ +ve parity Nj=2 Wilson 
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